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Abstract 



We study parametric resonance of interacting waves having the same wave 
vector and frequency. In addition to the well-known period-doubling instabil- 
ity we show that under certain conditions the instability is caused by a Hopf 
bifurcation leading to quasiperiodic traveling waves. It occurs, for example, 
if the group velocities of both waves have different signs and the damping 
is weak. The dynamics above the threshold is briefly discussed. Examples 
concerning ferromagnetic spin waves and surface waves of ferro fluids are dis- 
cussed. 
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Parametric resonance is an instability phenomenon that is responsible for the excitation 
of modes of a weakly damped system due to periodic modulation of parameters which 
determine the frequencies of the modes The prototypical example is the vertically 
driven pendulum where the modulated parameter is the acceleration of gravity ||. In 
spatially extended systems these modes are waves. Experimentally parametric resonance 
is well-known for spin waves (Suhl instabilities) f||-|5| and surface waves of fluids (Faraday 
instability) @-g. 

Let us briefly summarize the well-known behavior of the one-wave case where only one 
wave type is involved. The ground state (e.g., the flat surface of a fluid) becomes unstable 
if the driving (i.e., modulation) amplitude exceeds some threshold which is proportional to 
the damping constant. The instability is caused by the excitation of waves fulfilling the 
first-order parametric resonance condition 

m = f , (i) 

where f2(k) is the frequency of the wave with wave number k (i.e., the dispersion relation) 
and uj is the modulation frequency. At the instability threshold a period-doubling bifur- 
cation occurs because the wave solution has a period which is two times the period of the 
modulation. That is, in the stroboscopic map (i.e., the map which gives the state at time 
t = 2(m + 1)tt/u provided the state at t = 2rmr /uj is known) a period-two orbit bifurcates. 

How does this one-wave picture change if two types of waves with different dispersion 
relations fli and Q2 are involved? There is in principle no change as long as the manifolds 
defined by Qj(k) = uj/2 do not intersect. At such intersection lines two different waves with 
the same wave vector fulfill the parametric resonance condition. Interaction leads to the 
well-known effect of hybridization where the pure waves are replaced by mixtures of them. 
Although hybridization avoids the intersection of the above-mentioned manifolds they still 
come very close together with a minimum distance proportional to the strength of interac- 
tion. Below we will see that under certain very general conditions a direct Hopf bifurcation 
occurs. It leads to quasiperiodic motion right at the onset of parametric resonance. 
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The wave lengths of waves which are excited due to parametric resonance are usually 
large compared to the atomic scales. Their dynamics are described by macroscopic field 
equations like the Navier-Stokes equations for fluids or the Landau-Lifshitz equations for 
magnets. Formally the fields can be put together into a single vector field u(r, t), called 
order parameter. In the limit of weak damping universal equations of motion so-called 
amplitude equations can be derived |] . They describe the spatio-temporal evolution of the 
amplitudes of waves fulfilling the parametric resonance condition. The amplitude equations 
are derived from a multiple-scale perturbation theory where the damping constant is the 
smallness parameter. For the field u we make the ansatz: 



where A± and B± are the slowly varying amplitudes of the right and left waves, respectively. 
The corresponding wave solutions are u± and v±. In general they depend on y and z. The 
index ± corresponds to the different types of waves. The cross point of their dispersion 
relations Q±(k) are denoted by k and luq, i.e., fl±(k ) = u . It is assumed that the modu- 
lation amplitude h, the detuning uj/2 — Q(k), the amplitudes of the waves, and the coupling 
strength between different types of waves are of first order in the damping constant. In the 
one- wave case the amplitude equations for A± and B± are well-known [^jp|,|lO|. Coupling 
between + and — waves lead to additional linear cross terms. Because of the hamiltonian 
nature of these coupling u± and v± can always be chosen in such a way that the linear part 
of the amplitude equations reads: 



u 



(r,t) = (A + u + + A_u_)e i(fca; - LJ ' /2) + 

+ (5+v+ + R_v_)e <(te+wt/2) + c.c. 



(2) 



A 



± — 



-T±A± + iha±B± + iee ±ix A 



(3a) 



B± = -T* ± B± - iha ± A ± + iee Tix B : 



(3b) 



with 



T ± = 1± + i(Sl ± -u/2), 7± >0, 



(4) 
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where r y±, Q±, and a± are real numbers denoting for each wave type the damping constant, 
the dispersion relation, and the strength of modulation, respectively, e and \ are the strength 
and the phase of the coupling, respectively. These parameters are in general functions of 
the wave number k. The amplitude and frequency of modulation are denoted by h and lu, 
respectively. 

The ground state A± = B± = becomes unstable if a solution A± = a±exp(At), 
B± = b± exp(Ai) exists with Re A > 0. For increasing h, A will eventually cross the imaginary 
axis defining a bifurcation point. Since the coefficients of (|f) depend on k these points 
become functions of k called neutral curves. There are two different types of neutral curves 
possible: (i) Neutral curves hs(k) where A = (soft-mode or stationary instability) and 
(ii) neutral curves hn{k) where ImA = ujh ^ (hard-mode or oscillatory instability). In 
the stroboscopic map they correspond to period-doubling bifurcation and Hopf bifurcation, 
respectively. The neutral curves hs and hu are the real solutions of 

Eh A -Fh 2 + G = (5) 

with E, F, and G given by 

E s = («+a_) 2 , (6a) 
Fs = |a+r_| 2 + |a_r + | 2 — 2a + a_e 2 cos 2%, (6b) 
G s = |r + r_ + e 2 | 2 , (6c) 

and 

E H = 1+1 .{a 2 + -a 2 ) 2 , (7a) 
F H = 2 7+7 _[(« 2 h - a 2 _ )(|r + | 2 - |r_| 2 ) + 
+2( 7+ + 7_)( 7 _a 2 , + 7 + a 2 )] + 

+ ( 7+ + 7 _) 2 (a+ + c?_ + 2a+a„ cos2x)e 2 , (7b) 
GH = 7 + 7-[(|r + | 2 -|r_| 2 ) 2 + 

+4(7 + + 7 -)( 7 -|r + | 2 + 7 + |r_| 2 )] 
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+|( 7+ + 7-)(r + + r_) e 



(7c) 



respectively. The solution hn is relevant only if the Hopf frequency uoh given by 




7-|r + | 2 + 7+ |r-| 2 -(7^ + 7+ « 2 _)/ l | 



(8) 



7+ + 7- 



is real. 



In the noninteracting case (i.e., e = 0) the characteristic polynomial (Jg) factorizes into 
two second-order polynomials and we are back to the one-wave case. The instability is always 
caused by a real A crossing the imaginary axes and consequently Hh does not exist. The 
neutral curve reads h$ = |r±|/a±. The instability threshold h c (i.e., the absolute minimum) 
is given by h c = h± = 7±/a±. 

The two- wave case (i.e., e ^ 0) is qualitatively much richer. Outside the hybridization 
area the one-wave picture is still valid. We assume that the threshold is given by the + wave, 
i.e., h c = h + . How does the type and the threshold of the instability change if we sweep 
through a hybridization area? First it can be rigorously proved that the threshold is always 
higher than h + . This can be easily understood because the coupling of a parametrically 
excited wave with another wave opens an additional dissipation channel. Thus the driving 
amplitude has to be increased in order to excite a wave parametrically. In the hybridization 
region where the coupling is strongest this effect is most pronounced. In the limit of strong 
damping or weak coupling this effect is of the order e 2 . In this limit direct Hopf bifurcations 
are not possible. 

For the other limit (i.e., strong coupling or weak damping) we have to distinguish be- 
tween two cases depending on whether the slopes of the dispersion relations (i.e., the group 
velocities c± = dfl±/dk) have the same sign or not. When both signs are identical the neutral 
curve hs has always two relative minima located near the parametric resonance condition. 
The reason for that is that the hybridized dispersion relations are monotonic functions of the 
wave number. Sweeping uj through the hybridization region one minimum increases whereas 
the other one decreases [see Fig. |T](a)] . There will be a point in the middle near u = 2ujq 
where both minima are of equal height. This is a co-dimension two bifurcation point where 
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a competition between two period-doubling bifurcations for two different values of k takes 
place. Assuming no fc-dependence of the group velocities c± and the parameters 7±, a±, e, 
and x we S e t approximatively for the threshold in the co-dimension two point 

, /o \ |c-7+ + C+7-I / n v 

ft Sc (2a;o) = . (9) 



'cia^ + c|ai — 2c + c„a: + a_ cos 2% 
The error is of quadratic order in the damping constants. Note that the threshold enhance- 
ment hsc(2uJo) — h + can be of the same order as h + [see Fig. |l|(a)]. 

If the group velocities have different signs the hybridized dispersion relations are no 
longer monotonic. Thus an interval around uo = 2u> exists where the parametric resonance 
condition can not be fulfilled. Therefore the threshold hs is much larger than in the case of 
equal signs (see Fig. |IJ) because it can not be of the order of the damping constants like (0). 

In the limit of strong coupling or weak damping the instability due to a period-doubling 
bifurcation can be beaten by an oscillatory instability (see Fig. [TJ). Its neutral curve hu is 
very well approximated by 

ir+ + r„ 



(10) 



I a\ + a 2 _ + 2a + a_ cos 2% 
The Hopf frequency ujh is in leading order given by e which is half of the hybridization gap 
at k = ho- The threshold 

hnM = . = ? + + 7 ~ (11) 



'«! + «: + 2a + a- cos2x 
occurs at 

n + (k Hc ) + n4k Hc ) = Lu. (12) 

Eq. ([II] ) is a good approximation only near uj = 2uj . Going away from this point hfj c (u) 
increases [see Fig. |l](b)] . It eventually disappears because hu does not exists far away from 
the hybridization area. 

What are the conditions for hue < ^5c? Because of too many parameters we can not 
give a complete answer. We will consider only two important case where an easy answer 



is possible. In the first case the group velocities c+ and c_ have different signs and hue is 
almost always less than h Sc because hn c ~ an d h Sc ~ 0(7°). For equal signs we get 

an easy conditions if only one wave is driven parametrically (i.e., a_ = 0). Here (|9|) and 
([□]) simplify to hs c = (7+ + y~\c+ / C-\) / a + and hue = (7+ + 7-)/«+, respectively. Thus 
the Hopf bifurcation occurs first if |c_| < |c + |. 

In a preliminary study we have also investigated the bifurcating solutions by including 
usual third-order terms in the amplitude equations |l|J^,|T0[| . The bifurcation at hn is either 
supercritical or weakly subcritical. The bifurcating solution is similar to the drift solution 
which occurs in a secondary instability in the one- wave case (5,0 ■ For larger values of 
h this solution becomes oscillatory unstable. The corresponding Hopf bifurcation leading 
to a torus solution is either supercritical or weakly subcritical. Note that for the order 
parameter u this solution contains three different frequencies. Further increase of h leads 
via a period-doubling sequence to chaos. 

In which physical systems can we expect direct Hopf bifurcation due to parametric res- 
onance? First we distinguish two broad classes of physical systems where hybridization 
of waves on macroscopic scales takes place. The first class is characterized by at least 
two different types of waves caused by different physical mechanisms, e.g., spin waves and 
sound waves in ferromagnets. The phonon-magnon system has already been studied exper- 



imentally in ferromagnets by parallel-pumping |TTJ. Parallel pumping means an additional 
high-frequency magnetic field parallel to the static one. In parallel-pumping only the spin 
waves are driven. Unfortunately the usual period-doubling bifurcation comes first because 
the sound velocity is larger than the magnon group velocity. But the threshold enhancement 
has been observed |Tl]] . We predict a direct Hopf bifurcation if only sound waves are excited 



parametrically. 

The second class of systems with hybridization is obtained by squeezing a wave with an 
anisotropic and/or nonmonotonic dispersion relation into a geometry with small transverse 
extension. Examples are ferromagnetic films in saturating magnetic fields and ferromag- 
netic fluids in narrow channels with a magnetic field perpendicular to the surface. The 



corresponding dispersion relations are anisotropic Q and nonmonotonic | l2| , respectively. 



The different types of waves in such systems are caused by a discrete set of the transverse 
component kj~ of the wave vector. For a wave with a lateral component k the dispersion 



relation Qj is the bulk one where k is replaced by yk 2 + kf 2 . Crossings between fij+i and 
Qj occur when \kj~ +l — kj-\ is not too large, i.e., the extension of the system in the transverse 
direction is not too narrow. The strength of interaction is strongly influenced by the kind 
of boundary conditions. 

We have calculated the threshold for a parallel-pumped, ferromagnetic, insulating film 
in a uniform magnetic field oriented perpendicular to the film plane [fTi}| . Fig. |2] shows 
an excellent agreement between a fully numerical calculation and an analytical one within 
the framework of the amplitude equations (|3]). The coefficients of the linear terms of the 



amplitude equations have been already calculated in the literature [|T^,|TIJ . On the left hand 
side of Fig. |2| we see the neutral curves in the hybridization area of the modes j — and j = 2. 
The Hopf bifurcation has clearly the lowest threshold. But unfortunately a non-interacting 
wave (here the mode j = 1 on the right hand side) has a slightly lower threshold. Very 
recently Kostylev et al. jT7| have shown in a fully numerical calculation that in tangentially 



magnetized films the Hopf bifurcation can have the overall lowest threshold. 

A promising physical system for experimental verification of the direct Hopf bifurcation 
is the Faraday instability of a ferrofluid in a static field perpendicular to the surface in order 
to get a nonmonotonic dispersion relation. The fluid should be in an annular container with 
a width which is roughly as large as, or slightly larger, than the wavelength of the surface 
waves at the local maximum of the dispersion relation. At the cross points of dispersion 
relations of waves with different transverse wave numbers the group velocities have often 
different signs. Thus the Hopf bifurcation appears first if the damping is not to strong. 
Without any detailed calculation we can decided whether a non-interacting wave has still a 
lower threshold or not. 

In this letter the influence of hybridization on the parametric resonance of waves has 
been studied. We have shown that under very general conditions the usual period-doubling 
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bifurcation is replaced by a direct Hopf bifurcation leading to quasi-periodic traveling waves. 
The Hopf frequency is roughly given by half of the hybridization gap. 
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FIGURES 

FIG. 1. Thresholds and critical wave numbers as functions of the modulation frequency near 
a hybridization region. The parameters are e = 0.06, x = 1-5, a+ = 4, a_ = 0.5, 7+ = 0.03, 
7_ = 0.02, =a; + 0.2(A;-A;o), and (a) Q + = uj + 0-7(k - k ), (b) Q + = uj - 0.7{k - k ). Solid 
(dashed) lines indicate relative minima of hs (hn)- The absolute minimum is denoted by thick 
lines. 

FIG. 2. The neutral curves for a parallel-pumped insulating ferromagnetic film of thickness 
d. The magnetic field is uniform and perpendicular to the film plane. The values of the exchange 
length /, the static field H, the Landau-Lifshitz damping g, and the driving frequency u are 
0.05<i, 1.1 • AttMq, 0.01, and 0.42 • AttjMq, respectively, where Mq is the magnetization and 7 the 
gyromagnetic ratio. The solid and dashed lines denote the analytically found neutral curves hs and 
Hh-, respectively. The squares and triangles are numerical results based on a Galerkin expansion 
& 
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Fig. 2 — F.J. Elmer 
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